Dependence of Deep Inelastic Structure 
Functions on Quark Masses 



A.V. Kisselev, V.A. Petrov 

Institute for High Energy Physics 
142284 Protvino, Russia 



Abstract 

We argue that the difference between the structure functions cor- 
responding to deep inelastic scattering with and without heavy 
quarks in the current fragmentation region scales at high Q 2 and 
fixed (low) XBj- The lower bound on charm contribution to the 
total structure function, F|(Q 2 ,x), is calculated and compared 
with the recent data on F£(Q 2 ,x) from HI Collaboration. 



1 Introduction 



Quite often mass effects in high-energy collisions are considered as some not 
very spectacular corrections which finally die off. Nonetheless, it appears that 
in e + e~ annihilation even such overall characteristics as hadron multiplicities 
are quite sensitive to the value of masses of the primary qq pairs Q . 

Recent considerations have shown that calculations based on QCD agree 
well with the data at high enough energy || and that they yield an asymp- 
totically constant difference between multiplicities of hadrons induced by the 
primary quarks of different masses. 

In this paper we study a similar effect in a deeply inelastic process M, 
[|]. As a by-product, we estimate heavy quark contributions to the total 
structure function. 



2 Calculation of quark mass dependence 

Let us consider, for definiteness, deep inelastic scattering of the electron 
(muon) off the proton. The hadronic tensor (an imaginary part of the virtual 
photon-proton amplitude) is defined via the electromagnetic current J /t : 

W^(p,q) = \{^f J d 4 zexp(iqz) < p\[J„(z) , J v {0)]\p > , (I) 

where p is the momentum of the proton, p 2 = M 2 , and q is the momentum 
of virtual photon, q 2 = —Q 2 < 0. 

A symmetric part of has two Lorentz structures: 

W, v = {-g, v + q -f^ Fl (Q 2 lX ) + 1 (p„ - g„££) (p v - ft*f) F 2 (Q 2 ,x), 

(2) 

where the structure functions F\ and F2 depend on Q 2 and on the variable 

x = 91 (3) 

pq + yj(pq) 2 + Q 2 M 2 

In what follows we will analyse the structure function F 2 of deep inelastic 
scattering with open charm (beauty) production at small x. In this section 
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we consider the case of one single quark loop with mass m q and electric 
charge e q . A general case will be discussed in Section 3. 

At small x a leading contribution to F 2 comes from one photon-gluon 
fusion subprocess 0: 



W,u = I j^-^iq^^m^i^d^ik^'^i^p), (4) 



d A k 1 

where k is the momentum of the virtual gluon, k 2 < 0. The tensor 
denotes an imaginary two gluon irreducible part of the photon-gluon ampli- 
tude, while T Q/3 describes a distribution of the gluon inside the proton. A 
quantity d a p is a tensor part of the gluonic propagator. 
Let us choose an infinite momentum frame 



P» 



M 2 M- \ 

P+^L ,0,P- — ). (5) 
4P' AP ) V ' 



Then the gluon distribution T al3 has to be calculated in the axial gauge 
nA = with a gauge vector = (1, 0, 0, —1) ||. One can take, for instance, 



n 



% + xPn (6) 



with x defined by Eq. (| 
From Eq. (0) we get 



1 -f, 

x 



3Q 

-9iiv + V^Vv 



Gog) 2 + Q 2 M 2 



= F^ a) + F 2 {b) . (7) 



Two terms in the RHS of Eq. ([F]), F^ and F 2 \ correspond to two tensor 
projectors, and p^p u . 

Note that the structure function F L = F 2 — 2xFi is completely defined 
by the term p^p^ and, thus, proportional to F 2 . 

By definition, the gluon distribution T af3 can be rewritten in the form 

r ^ = J- E 5 (P + k ~ Pn) < P\I 9 M \n >< n|/|(0) \p >, (8) 

'±71 n 

where l 9 a is the conserved current. Both \p > and \n > are on shell states 
that result in 

k a T aP = 0. (9) 
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From an explicit form for C?£ (see Ref. Appendix I, for details) one 
can verify that it obeys the same condition: 

= o. (io) 

Equations @ and (|T0|) allow us to simplify expression (|]) and get (r = 

a, b): 

V = / (04^cS(ff,*;^)r^(*,P), (ii) 

with the notations 

~ (pg) 2 + Q2 M 2 ^ 
The tensor r Q/3 can be expanded in Lorentz structures 



knj \ knj k 2 \ knj \ knj k 

with r< = Ti(k 2 ,M 2 ,pk). 

Let us consider a contribution of the invariant function Ti into the struc- 
ture function F 2 ([O]). With the accounting for @ and ([[(]) we obtain 

I F M = P i) d ± Q f X) dl 2 l-l 2 x 2 /Q 2 z 2 (r) (Qp_ m\ x\ d 2 
x 2 ? i 2 7 / 2 1 + M 2 x 2 /Q 2 I / 2 ' Z 2 ' J dln/ 2 ^ 1 ' J ' 

(14) 

where 

I 2 = -k 2 > 0, (15) 

z = - (16) 
pn 

and 

= f^~ z (17) 
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Here we used the notation: 



C (r) = _ q aP C (r) 



a/3- 



(18) 



To be more correct, one has to write z > x(l+4m 2 /Q 2 ) and I 2 < Q(z/x) 
4m 2 z/(z — x) in flT4|), but me neglect power corrections 0(m 2 /Q 2 ). 
In Eq. (|14] ) the gluon distribution, g(l 2 ,z), is introduced: 



If we use the new variable 



2(2 



(19) 



pA; + V(pA;) 2 - fc 2 M 2 



(20) 



instead of kj_, we will arrive at the expression 



A thorough analysis shows, however, that the main contribution to F 2 at 



small x comes from T 2 and I\ in 
for details): 



and F 2 is given by the formula (see 



x 



AT. 



dz 

z 



\z/x) 



d 



J 2 ' I 2 ' z) <91n/ 2 

")2 >>-} 2 rr>\ ^ 



Q^rn^x 



/ 2 ' / 2 ' 



9 In/ 2 



G(/ 2 ,^) 



G(/ 2 ,z) 



.(22) 



As we are interested in a calculation of the difference of the structure 
functions corresponding to the massive and massless cases, we preserve those 
terms in C^ r ' which give a leading contribution to AF 2 . In || we have 
calculated the functions C^ a > in lowest order in the strong coupling a s : 



C {a) (u,v,y) 



Air 

a 

TV 



{[(1 - y) 2 + y 2 ]L(u, v, y) - [(1 - yf + y 2 - 2v]M(v, y)-l}, 



y{l-y)M{v,y), 



(23) 



where 

u(l - y) 



L(u,v,y) = In 



y[v + y(l - 



M M = f ;») ■ (24) 
As for the gluon distributions, they are given by the formulae: 

G = sk/^/f(f-*)(«' + J)[r.(' a .f) + r 1 (« a ,a P5) 



'V 2 f ../.. 2 r 2 ' 



\2 



4 £2 



-r 2 (/ /2 ,0 + r 4 r,0 



(26) 



The analogous expressions for the functions are the following [0: 
C^W.!/) = ^ y {2y[(l-2y)(l-y)-v)L(u,v,y) 

2,71 U 

o 

+ (1 - y)[{l - yf + y 2 - 2v]M(v, y)} + ^y(l - y), 

6V(u,v,y) = -^V(l-l/) S #,»)- (27) 
7r it 

It may be shown that the leading contribution to AF 2 comes from the 
region I 2 ~ m 2 , k 2 = —I 2 being the gluon virtuality. Then one can easily 
see from (p4j) and (p7|) that the first two terms in are suppressed by 
the factor k 2 /Q 2 with respect to C^ a \ while the third terms in do not 
contribute to the difference C^| m= o — C^lm^o- 

In the leading logarithmic approximation (LLA), only the function L 
remaines in Eqs. (|23|), which results in 



1 d „ , „ a q r dz 



xdhiQ 2 



F 2 (Q 2 ,x) = ^ [^pJ*)g(Q 2 ,z), (28) 

Z7T J Z \Z / 



where P qg {z) is the Altarelli-Parisi splitting function and G(Q 2 ,z) is the 
gluon distribution in LLA defined by Eq. (|25|). 
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It is clear from (g2|) that AF 2 = F 2 
tities (r = a, b) 



m= o — Fzlm^o is defined by the quan- 

(29) 



ACW («, v, y) = C« (u, 0, y) - C« (u, v,y). 
By using Eq. (|23| ) we obtain the important result 



while from Q27| ) we get 



AC (o) 
AC (a) 



AC< 6 ) 
AC< ft > 



AC^ )2 /), 



-A<7«(«,v), 
-AC^,y). 



(30) 



In this, we have 

So, we get M 
1 



AC (a) , AC (a) | 



-fc 2 - 



fc 2 



(31) 



(32) 



,r 



AF 2 (Q 2 ,m ,x)| Q2 



+oo e g 



1 T <i/ 2 



AC 



m 2 x\ 9 



Of! 



I 2 ' z I d\n P 



ml x\ d 



+ AC — r, - 



/ 2 z ) d\\il 2 



G{l 2 ,z) 



Here 



AC7(«,y) = ^[(l-^ + ^ln 



2/(1 -y) J u + y(l-2/) 



AC(^) 



a 



y(l — y) 

7T v + y(l-y) 



(34) 



with G(/ 2 , z) and G(Z 2 , 2) being defined by Eqs. (H) and ©. 

The integral in Z 2 (|33|) converges because of condition fl32|). Contributions 
from AC^ and AC^ are suppressed by the factors (m 2 /Q 2 ) InC; 2 and can 
thus be omitted. 
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Let us consider the gluon distribution G ( p^) . At small z the leading 
contribution to G(l 2 , z) comes from the region z <C £, and we have 

G{l 2 ,z) ~G{l 2 ,z). (35) 

Taking expression fl35|) into account, the structure function F 2 ( p2[ ) has the 
following form at low x (with the term of the order of k 2 /Q 2 and m 2 /Q 2 
subtracted) 

1 2 } dz Q { f X) dl 2 n ( Q 2 m 2 q x\ d a 

where 

C(u, v, y) = ^{[(1 - y) 2 + y 2 ]L( M , v, y) - [(1 - 3y) 2 - 3y 2 - 2v)M(v, y) - 1}. 

(37) 

As for the difference of the structure function, we obtain the following 
prediction 

(38) 

where 

AC^,y) = ^[(l-y) 2 +2/ 2 ]{ln 



1 + 



ufl — u). 



■(l-2y) s 



i> + v(l — v) 



■ (39) 



3 Relation between measurable structure func- 
tions 

Up to now, we considered those contributions to F 2 that came from the 
quark with electric charge e q and mass m q , F 2 \ m ^ . Then we have taken the 
analogous contributions from the massless quark with the same e q , F2|m=o, 
and have calculated the quantity -F 2 | m=0 — i^lm^o- 
The total structure function F 2 has the form 

F 2 (Q 2 ,x)=]Te 2 F|(Q 2 ,*), (40) 
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where the functions F 2 are introduced (g = u, d, s, c, b). 

The structure functions describing open charm and bottom production 
in DIS, F£ and F 2 respectively, are related to F 2 C and F 2 by the formulae 

4 ~ 

F b 2 = \F b 2 . (41) 

At low x one can put {m u = = m s = is assumed) 

i> = F 2 d = F° = F 2 (42) 

and define the difference between heavy and light flavour contributions to 
F 2 : 

AF 2 C = F 2 -F 2 C , 

AF 2 fe = F 2 - F*. (43) 

Notice that there are the functions F 2 and F 2 that have been calculated in 
the previous section (see Eqs ( |36l) and (|38|) ). 

From Eqs. fl38|) and (^3) one readily obtains that a linear combination 

£ a (Q 2 ,x) = F 2 {Q\x) + aF<(Q 2 ,x) - (4a + 11)F 2 "(Q 2 , x) (44) 

scales at Q 2 — > oo and arbitrary parameter a. In terms of AF 2 introduced 

lim £ a (Q 2 ,x) = -^(l + a)AF 2 (m 2 ,x) + ^(4a + 10)AF 2 (m2,x). (45) 

Q 2 ^oo 9 9 

Let us now represent the function F 2 (KB) in the following form 



f^Jd V C( V ,y)G'(Y-rj^y (46) 



where we denote 

r=in S (47) 
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and introduce the variable r\ = ln(k 2 /Ql). Here G' means the derivative of 
G(Q 2 ,x) with respect to the variable InQ 2 . 
Analogously, we get from (EH) 



1 



-AFf = / ?L 
x J y 



1 dy 



x 



dr l AC(r l ,y)G , [Y m - v ,- 

\ y 



(48) 



with 



m„ 

Y,„ In ' 



yQl 



(49) 



Here r] = ln(m 2 /k 2 y(l —y)) — ln(m 2 /k 2 y) (remember that we consider small 
x). 

The expression for AC is given by Eq. (^) and, in terms of the variables 
rj and y, looks like 



a. 



^C = ^[(l-y) 2 + y 2 } 



ln(l + e^ 



1 + e" 



(50) 



As for the expression for C, it has to be defined via relation (O) and 
exact formulae (HI) taken at m = 0. The result is of the form 



C(v,y) 



2^ 



(61) 



where 



U 
V 



'l-4y(l-y)e-n, 

(l-y)[y+(l-y)e-'>] (l 

It is clear from ( |5T^) that 

A% 2 /)>0 



(52) 



(53) 



for — oo < rj < oo, < y < 1 and AC(r),y) is negligible at 77 < (see 
Figs. la-Id). 

Moreover, the quantitative analysis shows that at most at y < 0.2, which 
is relevant for small under consideration, one has 



C(r],y) > AC{r),y), 77 > 0, 



(54) 
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(see Figs. 2a-2d). Neglecting the small contribution to F 2 from the region 
r] < and taking into account that dG(Q 2 , x)/d\nQ 2 > at small x (cf. 0), 
we thus conclude 

AF|(m 2 ,x) < F 2 (Q 2 ,x)\ Q 2 =m2q . (55) 

From Eqs. ([45]), ([55]) we obtain the following inequality which holds for 
-2.5 < a < -1 

0<S a (Q 2 ,a;)| Q2>>1GcV2 < + a)(F 2 - F 2 C - F 2 6 )(Q 2 , x)\ Q 2 =m , 

+ i(5 + 2a)(F 2 - F 2 C - F 2 6 )(Q 2 , x)| Q2=m2 (56) 

At the endpoints a = —2.5 and a = — 1 we get 

(f 2 -2.5F 2 c -F 2 6 )(Q 2 ,x) < (F 2 -F 2 c -F 2 b )(g 2 ,x)| Q2=m2 , 
(F 2 -F 2 C -7F 2 6 )(QV) < (F 2 -F 2 c -F 2 6 )(Q 2 ,x)| Q2=m2 . (57) 

Data on the total structure function F 2 for Q 2 between 1.5 GeV 2 and 
5000 GeV 2 and x between 3 x 10" 5 and 0.32 are now available f7j. As for 
the charm structure function, there are recent data on F 2 at Q 2 = 12 GeV 2 , 
25 GeV 2 and 45 GeV 2 with rather large errors 0. 

Using the first of the inequalities ( |57"|) we get (assuming F|(m 2 , x), F^im 2 , x) 

o(cf. !)) 

F 2 C (Q 2 , x) > 0.4 [F 2 (Q 2 , x) - F 2 6 (Q 2 , x) - F 2 (m 2 , x)] . (58) 

Taking use of the available data on F 2 (Q 2 ,x) 0] and neglecting F 2 (as 
F2/F2 reaches at most 2 3% at HERA) we estimate the lower bound ac- 
cording to (|58|) . The result is exhibited in Figs. 3a-3c with m 2 = 2.5 GeV 2 . 
Experimental data on F| at several values of Q 2 and x are taken from Ref. [H]. 

Tables 1-3 present the result of our calculations of the lower bounds on F 2 
at the same Q 2 and x. Three different values of F^ for each Q 2 , x correspond 
to m c = 1.3 GeV, m c = 1.5 GeV, m c = 1.7 GeV, respectively. 
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< X > 


F 2 C (theor.) 


F 2 C (exper.) |j 




.0008 


0.173 
0.161 
0.145 


0.211 ± 0.049 t °oil 


.0020 


0.137 
0.128 
0.116 


0.263 ± 0.036 ±g;gg 


.0032 


0.120 
0.112 
0.101 


0.190 ± 0.054 t° °oll 



Table 1. The lower bounds on F 2 C (Q 2 , x) for Q 2 = 12 GeV 2 . 



< X > 


F 2 C (theor.) 


F 2 C (exper.) § 


.0008 


0.258 
0.247 
0.231 


0.324 ± 0.099 ±g;gg| 


.0020 


0.205 
0.196 
0.184 


0.253 ± 0.069 ±g;gg 


.0032 


0.179 
0.172 
0.161 


0.222 ± 0.066 t° °oit 



Table 2. The lower bounds on F 2 c (Q 2 ,x) for Q 2 = 25 GeV 2 . 



< X > 


F 2 C (theor.) 


F 2 C (exper.) § 


.0020 


0.258 
0.249 
0.237 


0.156 ± 0.070 ±8;gl| 


.0032 


0.226 
0.218 
0.207 


0.275 ± 0.074 ±8:gta 


.0080 


0.165 
0.160 
0.152 


0.200 ± 0.064 
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Table 3. The lower bounds on F 2 C (Q 2 , x) for Q 2 = 45 GeV 2 . 

These estimates of F^ agree with the recent data on the charm contri- 
bution to F 2 ||. Our inequalities are also in agreement with the 
results of Ref. JT(| where the ratio F^j F 2 was estimated. For a detailed com- 
parison of our predictions with the data, an improved measurement of the 
charm component F 2 is required. 

Conclusions 

In this paper we have demonstrated that the lowest-order quark loop contri- 
butions to the structure functions at small x contain mass-dependent terms 
which scale at high Q 2 . This effect can be observed experimentally, and we 
predict theoretical bounds for the corresponding contributions from c-quarks 
(see Eqs. (|56| ) and (|5TD, Figs. 3a-3c, Tabs. 1-3). 
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Figure Captions 

Figs. la-Id: AC(r), y) as a function of the variable 77 at several fixed values 
of y. 

Figs. 2a-2d: C(r],y) (continuous curves) and AC(rj,y) (dashed curves) as 
functions of the 

variable 77 (77 > 0) at several fixed values of y. 

Figs. 3a-3c: The lower bounds on F^IQ 2 , x) (cotinuous curves) together 
with results 

from HI H (open and closed circles). 
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